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Four-Loop QCD Corrections to the p Parameter 
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The four-loop QCD corrections to the electroweak p parameter arising from top and bottom quark 
loops are computed. Specifically we evaluate the missing "non-singlet" piece. Using algebraic 
methods the amplitude is reduced to a set of around 50 new master integrals which are calculated 
with various analytical and numerical methods. Inclusion of the newly completed term halves 
the final value of the four-loop correction for the minimally renormalized top-quark mass. The 
predictions for the shift of the weak mixing angle and the W-boson mass is thus stabilized. 

PACS numbers: 14.65.Ha, 14.70.Fm, 12.38.Bx 

Electroweak precision measurements and calculations 
provide stringent and decisive tests of the quantum fluc- 
tuations predicted from quantum field theory. As a most 
notable example, the indirect determination of the top 
quark mass, m^, mainly through its contribution to the 
p parameter , coincides remarkably well with the mass 
measurement performed by the CDF and DO experiments 
at the TEVATRON 0. Along the same line, the bounds 
on the mass of the Higgs-boson depend critically on the 
knowledge of m* and the control of the top-mass depen- 
dent effects on precision observables. 

A large group of dominant radiative corrections can 
be absorbed in the shift of the p parameter from its low- 
est order value pBorn = 1- The result for the one- loop 
approximation 



5p = ixt 
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hence quadratic in rrit, was first evaluated in Q and 
used to establish a limit on the mass splitting within 
one fermion doublet. In order to make full use of the 
present experimental precision, this one-loop calculation 
has been improved by two-loop ^ H, 1^ and even three- 
loop QCD corrections Also important are two-loop 

H Ealll El [l3j and three-loop ^14, M electroweak ef- 
fects proportional to and x^, respectively, and the 
three- loop mixed corrections of order a^Xj |l5|. 

An important ingredient for the interpretation of these 
results in terms of top mass measurements performed 
at hadron colliders or at a future linear collider is the 
relation between the pole mass and the MS-mass def- 
initions, the former being useful for the determination 
of mt at colliders, the latter being employed in actual 
calculations and in short-distance considerations. To 
match the present three-loop precision of the p param- 
eter, this relation must be know in two-loop approxi- 



mation 0, 0, 0, and for the four- loop calcula- 
tion under discussion the corresponding three-loop result 
must be employed. 
For fixed pole mass of the top quark, the three-loop 
result leads to a shift of about 10 MeV in the mass 
of the VF-boson as discussed in 0i US- (This apphes 
both to the pure QCD corrections and the mixed QCD- 
electroweak one.) Conversely, the corresponding shift of 
the top quark pole mass amounts to 1.5 GeV. (Similar 
considerations apply to the effective weak mixing angle 
and other precision observables.) These values are com- 
parable to the experimental precision anticipated for top- 
and H^-mass measurements at the International Linear 
Collider |3| • In addition, there exists a disagreement (on 
the level of 3(t) between the values of the so-called on- 
shell week mixing angle, sin^ Ow, as measured by NuTeV 
collaboration in deep-ienalstic neutrino scattering |25| 
and as obtained from the global fit of the Standard 
Model to the electroweak precision data js^l • 

From all these considerations an improvement of the 
theoretical accuracy seems, therefore, desirable. This, 
however, requires the evaluation of four-loop QCD cor- 
rections to the p parameter, the topic of the present work. 
The shift in the p parameter is given by 
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where 11^^'^ (0) are the transversal parts of the W- and 
Z-boson self-energies at zero momentum transfer, respec- 
tively. The calculation is thus reduced to the evaluation 
of vacuum (tadpole) diagrams. 

The IF-self-energy receives contributions from the cor- 
relator of the "non-diagonal" i-6-current only. Contribu- 
tions to the Z self-energy originate only from the "diago- 
nal" axial current correlator induced by top quark loops. 
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The vector current part vanishes due to current conserva- 
tion, the bottom quark is taken as massless. The nonva- 
nishing parts of n^(0) are conveniently decomposed into 
non-singlet and singlet pieces characterized by Feynman 
diagrams where the external current couples to the same 
and to two different closed fermion lines, respectively. In 
three-loop approximation the singlet-piece is larger than 
the non-singlet piece by nearly a factor twenty. This has 
motivated the authors of to evaluate, in a first step, 
the four-loop singlet piece. The strategy employed in that 
paper was based on the algebraic reduction of the ampli- 
tudes to a small set of master-integrals with the help of 
the Laporta algorithm [isl Efll l an approach which was 
recently also applied to the evaluation of the two lowest 
non- vanishing Taylor coefficients of the vacuum polariza- 
tion l3Ct I3II and to the decoupling of heavy quarks in 
QCD [sdlSSl ls^. both in four-loop approximation. 

In order to complete the evaluation of the four-loop 
QCD corrections to the p parameter, 11^(0) and the 
non-singlet parts of Ilf (0) are required. The evaluation 
of nf (0) is fairly straightforward: the input diagrams 
have been generated with QGRAF [s^; for the algebraic 
reduction to master integrals an efficient program has 
been constructed Hi which relies on F0RM3 and FER- 



MAT mUm^. Furthermore, the full set of the corre- 
sponding master integrals is available with high precision 

The evaluation of 11^(0), however, requires the knowl- 
edge of a sizeable number of new master-integrals, a ma- 
jor part of them (around 40) nontrivial to evaluate pre- 
cisely. The master integrals can be chosen in many differ- 
ent ways. As discussed in ji^l the choice of a so called "e- 
finite basis" leads to integrals particularly suited for the 
evaluation through Fade approximations. On the other 
hand, topologies with eight lines or less are conveniently 
calculated through difference equations. In the present 
paper we therefore employ a combined approach, which 
makes use of difference equations [i^ E3| to evaluate the 
simpler topologies, i. e. those with up to eight lines, and 
a semi-numerical method based on Fade approximations 
El El El El El- There, a suitably chosen line of the 
four-loop vacuum diagram is cut, the large- and the small 
behaviour of the resulting three-loop propagator are 
calculated analytically El El El Ell , the function in 
the whole region is represented by Fade approximations 
and the remaining integration is performed numeri- 
cally (see Fig. 1). 
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dq F{q^) 



Fig 1. Symbolic description of the Fade method: One line of the vacuum integral is cut, the resulting propagator 
is represented by a Fade approximation and integrated num erically. 



An estimate of the numerical uncertainty is obtained 
by comparing different Fade approximations based on 
the same input information from the large and small g^- 
region, or by increasing the input information through 
inclusion of more terms from the high and the low q^ 
region. Furthermore, in all cases at least two different 
lines were cut to check the consistency of the results. A 
detailed discussion of the various applications, character- 
istic examples and comparisons with analytic results, e.g. 
for the lowest moment of the polarization function can be 
found in Ell- 
in contrast to the applications discussed in earlier pub- 
lications [il El Ell I massless cuts unavoidably arise in 
some of the relevant diagrams and a generalization of 
the method is required: In addition to the introduction 
of a suitably chosen function needed for the subtraction 
of the high energy logarithms, another function is em- 
ployed for the subtraction of the logarithms arising from 
the massless cut in the low energy limit. 

The result for the shift in the p parameter can be cast 



into the following form: 

<5p^=3x.^(^)^5pP (3) 

i=0 

Here Xt is expressed in terms of the MS quark mass 
rrit = at scale /i = nit, and as, defined in the MS- 

scheme for six flavors, is chosen at the same scale. The 
normalization factors are such that Sp^^ = 1. 

For the four- loop non-singlet result, decomposed ac- 
cording to the various color structures and the n/- 
dependence, we find: 

5p|^(non-singlet) = (4) 
1.5211 C|. + 1.2363 C|Ca 
-2.3132 ClTni- 4.5962 C| T nu 
-F0.7438 CfCI - LHQbCpCATni 
-h2.5037 Cf Ca T rih + 0.4681 Cp nf 
+0.6880CFT'^nl+ 0.8495 Cf n^n;. 
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with Cf = (N^ - l)/(2iVc), Ca = Nc, and T = 1/2, 
where Nc = 3 is the number of colors, nf denotes the 
number of active (hght plus heavy) quark fields, with 
Uf = ni + n}i. This result has been also independently 
obtained with the help of direct application of the Pade 
approximation method like it was described in |45L 
for the lowest Taylor coefficients of the vacuum polar- 
ization. We find agreement for all color structures with 
the relative accuracy varying between 0.4% and 4%, and 
confirm the result for the singlet contribution. 

Setting n/i = 1, n/ = 5 and the color coefficients to 
their natural values, we find 

Spf^ = (5pfiS(gij^gig^) _^ ,5p|^(non-singlet) 

= -3.2866+ 1.6067 = -1.6799, (5) 

where we have also displayed the result of 27] for the 
singlet piece. The singlet piece is still larger than the 
non-singlet piece by a factor two. Nevertheless, the hi- 
erarchy is less pronounced than in the three-loop case. 
Numerically, the overall correction looks small, just as in 
the two- and three-loop case. However, if the result is 
expressed in terms of the pole mass, a major shift origi- 
nates from the large correction in the pole-MS relation: 

SpP"'' = -93.1501. (6) 

For fixed top mass, this corresponds to a shift of around 
2 MeV in the M^-boson mass, well below the precision 
anticipated for future experiments. 

In conclusion, the full OiXtO?^^ contribution to the p 
parameter proves to be small and the result based on the 
three-loop calculation is stabilized. 
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Note added. 

The results of our calculations have been recently con- 
firmed in the independent work [51i |. 
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